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TECHNICAL NOTE 


GRAPHICAL SOLUTION OF SOME AUTOMATIC -CONTROL PROBLEMS 
INVOLVING SATURATION EFFECTS WITH APPLICATION 
TO YAW DAMPERS FOR AIRCRAFT 
By William H. Phillips 

SUMMARY 


A graphical method is presented for determining the motion of a 
freely oscillating system of one degree of freedom stabilized by a 
controlling device which applies control force in proportion to the 
displacement of the system, to its rate of change of displacement, or 
both. The contralling member is assumed to have limitations on its 
maximum deflection and on its maximum rate of movement. Several 
examples are presented to illustrate the method. 

From these examples, it is sho'wn that, at sufficiently small 
amplitudes, the period and dajtigjlng of the system correspond to those 
provided by linear operation of the control whereas, at very large 
amplitudes, the period and damping approach those of the uncontrolled 
system. At Intermediate amplitudes, if the deflection of the control 
is limited, a smooth transition between these two conditions of period 
and daBQ)ing takes place whereas, if the rate of control movement is 
limited, the damping may be reduced below that of the uncontrolled 
system. In some cases, limiting the rate of control movement may 
produce instability over a range of an^litudes. 

If the control produces primarily an Increase in damping, the 
control remains effective in producing damping even at amplitudes 
several times that at which saturation effects are first encountered. 
This effect may be oiseful in reducing the power requirements of yaw 
dampers for aiiplanes . 


INTRODUCTION 


Autopilots and other automatic control devices are frequently 
designed on the basis of linearized theory. In practice, however, the 
limits of lineaj operation of these devices may be frequently exceeded. 
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In some cases, "tlie use of a system designed to exceed its linear range 
may te advantageous. For exanple, recent flight tests of devices 
intended to in 5 ) 2 rove the damping of the lateral oscillations of 
airplanes have shown that large Increases in damping are required only 
at small amplitudes of oscillation. For this reason, it may he 
desirable to reduce the power and force output of the controlling 
device by limiting either its deflection range, its maximum rate of 
movement, or both, in order that the device may be made as small and 
as 1 •! gb-h as possible. In practice, such a control would frequently 
be required to operate in the range of amplitudes where these limitations 
are reached. Calculation of the effects of the device under these 
conditions is therefore of interest. 

In electrical work, effects resulting from operating beyond the 
linear range of a device have been called saturation effects. An 
analysis of the motion of a control system when saturation effects 
are involved may frequently be accomplished by conventional methods 
but the process is tedious, particularly when a number of initial 
conditions are to be investigated. Problems of this type may be 
conveniently solved by Tise of an electric analog computer (refs. 1 
and 2 ) provided such equipment is available. If an analytical solution 
is required, or if an analytical check on the results obtained from an 
analog comptcter is desired, a simple graphical procedure for solving 
such problems may be of interest. 

A method frequently found convenient by previous investigators 
for analyzing nonlinear systems is the use of the phase plane. This 
plane, which consists of a plot of the velocity of the system against 
its displacement, has been tised by K. Minorsky (ref. 3) others for 
the study of problems involving nonlinear force variations. The same 
method has been applied by Herbert K. Weiss (refs. 4 and 5) to a study 
of relay servomechanisms which involve discontinuous force variations . 

By a slight modification of the usual phase-plane technique, Irmgard 
Fliigge-Lotz (ref. 6) has obtained solutions for the motion of an oscil- 
latory system such as an aircraft controlled by discontinuous or on-off 
controls. In the present report, methods somewhat similar to those used 
in reference 6 are applied to the solution of some problems Involving 
saturation effects. 

The following analysis is applied to the case of a freely 
oscillating system with one degree of freedom stabilized by a controlling 
device which applies control force in proportion to the displacement of 
the system, to its rate of change of displacement, or both. The 
controlling member is assumed to have 3jjnltations on its maxlmijm 
deflection and on its TnnxiTTniTti rate of movement. Graphical solutions 
for several problems of this type are presented to illustrate the 
method . 
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SYMBOLS 


A, B 

C 

D 

Iz 

N 

R 

r 

s 

As 


t 

X, y 

Sr 

e 

A 0 

? 


constants (see eq. (9)) 


constant 




B^ 


differential operator, d/ds 
differential operator, d/ds^^ 
moment of inertia in yaw- 
yawing moment, or ratio 

maximum rate of rudder n£»vement in terms of s 
radius (in polar* coordinates) 
nondlmensional measure of time, o^t 
change in nondlmensional time 

nondlmensional meas-ure of time for linear range of operation 
of control, sl/1 + N 


■> "(A 


^ 5 r 




time 


recteuigular coordinates (see eqs. (17) and (18)) 
rudder angle 

angle ( in polar coordinates ) 
change in angle 0 

damping ratio of uncontrolled airplane 

damping ratio for linear range of operation of control, 

86 

25 + N ^ 

8Dilr 


2 \/ 

86 ^ 
1 + N — ^ 

V 

8 ilr 




angle of yaw 
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0 ) 

‘“c 


N, - ^ 

Sn 


N6r = 


phase £ingle, tan“^ ^ 

£ 


damped natural frequency of xincontrolled aiiplane in terms 


of time unit s, \/l - ^ 

dan 5 >ed natural frequency for linear range of operation of 


control in terms of time unit s^; 






\indamped natural frequency of uncontrolled airplane^ 
radians /sec 

undamped natural frequency for linear range of operation of 

^sec 


control, 05^^ 1 + N , radians/s 


B5t 


Sijr 


Subscripts; 

0 initial value or value at start of an Interval 

max maximum value 

A dot over quantity denotes differentiation with respect to 
actual time. 


ANALYSIS 


Theoretical development.- An extensive analysis of the motion of 
TniRRilaa nndRr the influence of discontinuous automatic controls is 
presented in reference 6. This reference enploys a modified phase- 
plane method and contains a detailed development of the theory Involved. 
The theoretical background required for the present analysis is very 
similar to that utilized in reference 6. In particular, the methods 
of reference 6 could be enployed in the region where the control is 
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either agednst its stop or moving at Its maxirnnin rate. The methods of 
reference 6 , however, are not convenient for handling the transition 
from the nonlinear to the linear range of operation. Somewhat different 
graphical procedures are employed herein to handle this problem. In 
order to derive these methods, a review of the theory starting with the 
basic equations of motion appears desirable, although some repetition 
of the initial development of reference 6 is involved. 

In analyses of nonlinear systems by the phase-plane method, the 
velocity of the system is ordinarily plotted against its displacement. 
The curve traced out by the system during its motion is called a 
trajectory in the phase plane. For a single-degree -of- freedom linear 
system with damping performing a free oscillation, the trajectory may 
be shown to be a spiral converging on the origin. This spiral is not 
a true logarithmic spiral but is somewhat distorted. It is shown in 
reference 6, however, that, if the values of velocity and displacement 
are plotted on a skewed system of axes, a true logarithmic spiral 
results . The use of a logarithmic spiral is very advantageous for 
graphical analysis. Inasmuch as all such spirals for a given value of 
the damping ratio are congruent. As a result, only one spiral for a 
given value of the damping ratio needs to be plotted. This spiral may 
then be made to pass through any given point in the plane to fit the 
desired initial conditions sin^jly by rotation about the origin. In 
o2rder to obtain the relations required for the graphical construction 
described herein, expressions are required for the displacement and 
velocity of the controlled system for various conditions of the control. 
These expressions are now derived. 

For purposes of illustration, the equations of motion are set up 
for an aiiplane oscillating in yaw and controlled by the rudder. The 
same analysis may be applied, however, to any single-degree-of- freedom 
system obeying similar equations of motion and laws of control. The 
equation of motion is as follows: 


Iz'i - 


In problems of this type, the number of parameters may be reduced 
by expressing the oscillatory characteristics in terms of a natural 
frequency and a daiiping ratio ^ which gives the fraction of 

critical daniping (ref. 7). These quantities for the uncontrolled 
system are defined as follows: 
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Also, let N = 
motion "becomes 


% 


In terms of these parameters, the eguation of 

(1) 


A farther feduction in the niimher of parameters may he obtained 
hy use of the nondimens ional measure of time s = cs^t* Let 

D = — j equation ( l) then becomes 
ds 

+ 2^Dt + t = -5rN (2) 


When the system is in its linear range, the rudder angle is 
proportional either to the angle of yaw, to the yawing velocity, or to 
a combination of these quantities. The law of control for the rudder 
in the linear range of operation is therefore 






(3) 


When the limits of ILnear operation are exceeded, the rudder angle is 
either constant (against a stop) or changing at a constant rate. The 
rudder angle in this case may be expressed by the equation 

^r = ^0 ^ 

where R represents the constant rate of rudder movement in terms of 
the nondimens ional time unit s . The solutions of the equation of 
motion for these ca^es are now discussed. 


First, the solution for the lin ear range of operation is considered. 
When equation ( 3 ) is substituted into equation ( 2 ), the following result 
is obtained: 




( 5 ) 


Mate a second change in the time variable. 
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= 8 + n! 




( 6 ) 


and let D„ = — — . Equation (5) then changes back to the form of 


ds 


equation (2), 


Dc V + 2^c®c^ + Tjr = 0 


where 


^c = 


"5 


, “Ss" 




(7) 


(8) 


The quantity represents the damping ratio of the controlled 

airplane. The natural frequency of the controlled airplane in radians 
per second is obtained from equations (l), (2), and (5) as 

I 

X * V 

Before writing down the solution of equation ( 7 ), the solution of 
equation (2) is presented for the case in which the rudder is against 
a stop or moving at a constant rate. Under these conditions, the 
right-hand side of the equation is not zero. The solution of 
equation (7), the right-hand side of which is zero, is then obtained 
from this result as a special case. "When the expression for the rudder 
angle (eq. (4)) is substituted into equation (2), the solution for the 
variation of t with s for specified initial conditions tI^q, (D^f)Q, 
and by, is as follows: 

ijr = - 2^E^ - BNs + e”^^(A cos cos + B sin css) (9) 

where 

A = 4ro + N(b^Q - 2 Cr) 

= \/l - 


03 
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The value of inlr may he expressed as a function of s as follows? 

Dilr = -EN + e~^®C-cuA sin cds + oiB cos (jus ) - 

+ RN^ (lO) 

The solution of equation ( 7 ), which applies when the control is 
operating in its linear range, may he obtained as a simplified case 
from expressions C9) and (lO), because in this case the form of the 
differential equation is the same hut the right-hand side of the 
equation is zero. The solution of equation (7) is obtained by setting 
and R equal to zero and by substituting and for ? 

and D in equations ( 9 ) and (lO). 


The derivation of the relations required for the graphical 
construction, in which the motion is represented by a true logarithmic 
spiral rather than by a distorted spiral, is now presented. From 
equations (2) and (4) 

+ 2^1>lf = -Tjr - ^8ro + Rs)n (ll) 

Substituting the value of ijr from equation (9) into equation (ll) 
yields the relation 

-D^ - 2^(D'1<’ + RN) = e"^®(A cos cds + B sin ojs) (l2) 


By rearranging equation (lO), the following ecjuation is obtained: 

D\}r + RN + ^rii’+ N^8 j.^ - 2^R^ + RNsj = 0 )e~^®( -A sin cos + B cos cjus ) (l3) 


In equations (l2) and (13), set 


where C = 

- 2^(Dt + EIO] 


A = C sin 0 
B = C cos ^ 

and = tan”^ Ecpiation (12) then becomes 

_f g 

= ooCe ^ (sin ^ cos cce + cos ^ sin css) 

= cuCe“^® sin (ous + 0) (l4) 
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Equation (13) becomes 

Dt + EN + ” 03Ce"^®(sin 0 sin os + cos 0 cos ods) 

= 03Ce“^® cos (toe + 0) (l5) 


The values on the right-hand sides of equations (lifr) and (l5) may be 
identified as the rectangular coordinates of a logarithmic spiral 
vhich, in polar coordinates, has the formula 


where 9 = ci£ + 0 and tq 


r = rQe 


-Ce/c 


CO 


(16) 


and 


In a set of rectangular coordinates, therefore, values of 
y eire plotted such that 


X 




-D^ - 2^(Dilr + mJ 


y 


= Dilr + RN 




X 


(17) 

(18) 


The resulting curve, which represents a trajectory of the motion, is a 
logarithmic spiral centered at the origin. In practice, the initial 
angle 0, which is related to the initial conditions of the problem, 
never hsis to be determined becaiise,. as mentioned previously, a 
logarithmic spiral may be rotated about the origin to pass through any 
point without changing its shape. A relationship given with equation (16) 
shows that the nondimenslonal time elapsing for any change in the angular 
coordinate 9 is given by 


The convergence of the spiral is determined by the exponential factor 
^/(o. ^Actually, only one variable is Involved because, as stated in 

connection with eq. (9), cu = 

The coordinate system for the logarithmic spiral, of course, 
differs from the velocity-displacement coordinate system ordinarily 
used in phase-plane studies . From the rectangular coordinates x and 
y, however, the values of Dt, and '>1'' may be determined as 

follows: First, eliminate ilf between equations (ll) and (18). The 

result is 


y = Dilr(l - 2?^) - + EN( 1 - 2^^) 


(19) 
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Equations (iT) and (19) express x and y in terms of D^l'' and D%. 
If the value of iWf is eliminated "between equations (l?) and (l9)^ 


_ D^\lf _ ( 1 - 2^^)x 
" 2^ 2^0i 


( 20 ) 


Equation (20) allows lines of constant yawing acceleration to he 

drawn in the x,y plane. These lines are straight lines with a slope 


1 - 2 ^ 
2^(0 * 


The line for = 0 passes through the origin. 


If the value of D ijr is eliminated between equations (l?) and (19), 


y 


= Dt + RN + 


CD 


( 21 ) 


Equation (2l) allows lines of constant yawing velocity D'4'' to be drawn 
in the x,y plane. These lines are straight lines with the slope C/cd. 
The line for Dilf = -RN passes through the origin. 

In order to determine the angle of yaw ilf, substitute 
from equation (U) into equation (17)- The result is an alternate 
definition of x, 

X .a>p*N(b^0 - 2 £r^ + RN^ (22) 


This equation shows that lines of constant are parallel to the 
y-axls. The grid of lines of constant 'if is moving across the x,y 
plane at a constant rate for the case in which the rudder is moving at 
a constant rate. This result is shown by the term RNs, which represents 
a shift in the origin of the grid of ilr lines by an amount proportional 
to the nondimens lonal time s . 


Beca\ise such a moving grid is inconvenient for graphical work, the 
value of ■'I'' may be obtained more easily for the case in which the 
rudder is moving at a constant rate by reading values of and Ir'^ 

from their fixed grids and using equation (ll) to determine the value 
of ■'1'^. For the case where the rudder angle is constant, however, the 
value of R is zero, and equation (22) reduces to 



(23) 


In this case, lines of constant ilr eire fixed lines , peirallel to the 
y-axls. The line for ijr = -R^rQ passes through the origin. 
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The geometrical relationship among the three sets of grid lines 
for reading Dilr, and ijr is shown in figure 1. This relationship 

may he derived from equations ( 20 ) , (2l), and (23). The grid lines for 


Di)f constant are rotated from the horizontal by the angle tan~^ the 

grid lines for D 'I'' constant are rotated from the vertical hy twice 

-1 2 ^ 


this angle, which equals tan* 


The spacing between grid 


1 - 

lines for D^''^ = 0 and D^ijf = 1 is shown to be one unit when measured 
along a line of constant Dilr. Similarly, the spacing between the lines 
for Dilf = 0 and Dijr = 1 is one unit when measured along a line of 
constant The spacing of the grid lines for reading ijr is such 

that these grid lines would paiss through the intersections of a grid 
line for Dilr constant and the various grid lines for These 

geometric relationships may be useful in constructing the grids for 
graphical work. The grids drawn in figure 1 are shown centered at the 
origin of the x and y coordinates, which is the origin of the spiral. 
In use, however, the line Dilr = 0 often does not pass through the 
origin of the spiral because, as shown by equation (21), the line for 
DiJr = -RN passes through the origin when the rudder is moving at a 
constant rate. The spacing and direction of the grid lines, however, 
are the same as those shown in figure 1. Similarly, the grid for 
reading ijr is often moved laterally with respect to the other grids 
while its direction is maintained paralled to the y-axis. 


Application of theory .- The results just presented are now 
summarized briefly and an attempt is made to show their physical 
meaning. The analysis considers the motion of a linear system of one 
degree of freedom iinder the influence of a linearly increasing force 
(corresponding to a constant rate of rudder deflection), a constant 
force (corresponding to a fixed rudder deflection), or zero force 
(corresponding to motion in the linear range of operation). In each 
case, the motion is represented by a logarithmic spiral in a certain 
system of coordinates designated the x,y plane. Values of lAlf, Di]f, 
and ilr may be read from suitable grids of lines superimposed on this 
plane. 


In the case of a linearly increasing force, the displacement of a 
single-degree-of-freedom system is known to Increase linearly with time 
after the transient oscillation has damped. The velocity of the system 
approaches a constant value, and the acceleration approaches zero. 

This behavior agrees with that given by equations (22), (2l), and (20), 
which show that the grid of lines of angle of yaw move across the x,y 
plane at a constant rate, whereas the spiral representing the motion 
converges to a certain value of yawing velocity and to a value of yawing 
acceleration of zero. 



12 


NACA TN 303 ^ 


Ttie sirapler cases of a constant force or of zero force may be 
shown to correspond to equations (22), (2l), and (20) in a similar 
manner. In the case of a constant force, the displacement approaches 
a constant value, and the velocity and acceleration both approach zero. 
In the case of zero force, all three quantities approach zero. The 
origins of the grids for reading Di|r, and ilf for this condition 

all coincide with the origin of the spiral. The tilt of the grid lines 
for reading values of Dijf and results from the fact that, when 

a slngle-degree-of- freedom system with damping performs a free oscilla- 
tion, the velocity leads the displacement by 90° plus the angle 

tan~^ ^ and the acceleration leads the velocity by 90° plus the same 

angle. These relationships were pointed out by R. K. Mueller in 
reference 8. The present analysis shows that this same concept can be 
extended to the case of the system xmder the influence of a constant 
force or a linearly increasing force, but that a special coordinate 
system is required to take into account the constant and linearly 
increasing terms in the solution for the transient motion as well as 
the term representing the damped oscillation. 

In order to make a graphical solution of the motion of the system 
for specified values of the variables, suitable spirals and grids 
should be constructed. One spiral is required for the linear range of 
operation, based on the damping ratio of the controlled system. 

Another spiral is req\iired for the saturated range of operation, based 
on the damping ratio ^ of the uncontrolled system. A suitable set 
of grids for each of the spirals should be prepared on separate sheets 
of tracing paper. These grids may then be superimposed on the spirals 
with the origins located to represent the various conditions of constant 
rudder rate, constant rudder deflection, or linear operation of the 
control, lie procedure for obtaining a solution consists in tracing 
the trajectory representing the motion of the system, with suitable 
shifts in the grids or changes in the spiral at points where the mode 
of operation changes from one regime to another. (These points are 
called "switch points" in ref. 6.) The correct procedure 'for making 
these shifts is best described by means of specific exan^les . Several 
such examples eire now presented. 


EXAMPLES 


Five examples are presented to illustrate different features of 
the method. These examples are summarized in the following table: 



NACA TCI 3054 


15 


Example 

Control type: Rudder 

angle proportional to - 

Limitation 

I 

Yawing velocity 

Rudder rate limited 

II 

Yawing velocity 

Rudder deflection limited 

III 

Angle of yaw 

Rudder rate limited 

IV 

Angle of yaw 

Rudder deflection limited 

V 

Angle of yaw and 
yawing velocity 

Both rudder rate and 
deflection limited 


Specific values for the parameters and Initial conditions are assumed 
for each example. No attempt is made to investigate a range of these 
parameters or initial conditions. 

Example I.- The following parameters are assumed: 



^ = 0.05 Cco = 0.9987) 
N = 0.3 
R = iO.125 

from which may he derived, from equation (8), 

= 0.2 (coc = 0.9798) 

The following initial conditions are assumed: 


% = -0.531 

(Dt)o = 0 

6 _ = -0.25 

from which may he derived, from equation (ll), 

(l^t)Q = 0.606 


T"vfo spirals should he constructed, corresponding to values of damping 
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ratio of O.05 and 0.2. A typical spiral, with the associated grids for 
reading D^ilr and D^lr, is shown in figure 2. A convenient and suffi- 
ciently accurate method of constructing a logarithmic spiral as a 
succession of circular arcs is presented in reference 6 (pt. 1, 
appendix I, p. 4l). 


In figure 2, the x- and y-scales are drawn on the same graph as 
the spiral. This representation is correct, inasmuch as the x- and 
y-scales are in all cases considered to be centered at the origin of 
the spiral. When the spiral is used, however, it must ordin^lly be 
rotated about its origin to pass through the desired initial point. 

If the spiral is rotated, the x- and y-scales do not rotate but remeiin 
fixed in direction. In practice, therefore, the values of x and y 
are read with a ruler rather than from scales marked on the same sheet 
as the spiral. 


The graphical work for example I is shown in figure 3- In this 
figure and succeeding figures in which graphical solutions are shown, 
the scales labeled x and y are drawn, for convenience, with their 
origin at the origin of the grid for and These scales should 

be considered to show only the direction and spacing but not the true 
position of the x- and y-scales, because the x- and y-scales are in all 
cases considered to be centered at the origin of the spiral, which is 
shifted in various parts of the solution. These scales therefore move 
when the origin of the spiral changes . 


In the linear range of operation, 6^. = 


The quantity 


^5 

- (called the "control function" in ref. 6) determines the ludder 
oDilr 

position called for by the controlling mechanism. Since initially the 


control function = 0 (i) =0 whereas = -0.25, the rudder is 

not in the position called for by the controlling mechanism and will 
run at its maxlinum rate in the direction to reduce the error between 
its existing position eind the position called for by the controlling 
mechanism. These conditions establish that the initial interval of 
the motion takes place with the rudder moving at a constant rate in 
the positive direction. The spiral corresponding to ^ = O.O5 should 
therefore be used. From formulas (20) and (2l), the origin of the 
spiral should be located at the intersection of the grid lines = 0 

and Dilf = -EN = -(0.125)(0.3) = -0.0375* The spiral may now be rotated 
so that it passes through the initial point (l>t)o = (E^'*1'')q = 0,6o6. 

The section of the spiral corresponding to the first interval may then 
be traced in as shown in figure 3(a)* The ends of the intervals, which 
correspond to changes in the direction of motion of the rudder or to 
changes in the mode of operation, are shown in figure 3 by circled 
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points which are nxnnhered to correspond to the number of the integral. 

The radial line from the origin of the spiral passing through the 
Initial point corresponds to s = 0. Values of D^ijr and Dijr may he 
read from the grids at successive values of the radial angle 0 measured 
from this radial line. For each value of 9, the corresponding value 

of nondimens ional time is determined from the formula s = — . If 

(a 

desired, the values of ijr may he determined at the same points hy 
substituting the values of Ihjf, and s into equation ( ll) . 

As the motion continues, the error between the rudder angle 
existing and that called for hy the control is eventually reduced to 
zero. At this point the Interval ends. This condition exists when 


8-_ = + Es 

^ ^0 



Make the substitution s 


Q 

— and solve for Ihlf. 
m ^ 

5r^ + 

D\lr = -2- 




Then, 


The curve showing Dijr as a function of 0 satisfying this relation 
is shown in figure 3 as a short dashed line. The intersection of this 
curve and the spiral determines the time at which the rudder reaches 
the angle called for hy the control. Substitution of a few values of 
0 into this equation will quickly show the values of 0 that will 
give values of Dilr close to those required to intersect the spiral. 

As a result, only two or three points of the curve in the vicinity of 
the intersection need to he plotted. 


An alternate method for deter minin g the time at which the rudder 
reaches the angle called for hy the control is to carry along an 
a\rxiliary plot of the rudder angle 6p and the control function 
Shr- 

Ihlr •:—-=• as functions of s. The intersection of these curves then 

determines the desired time. This type of plot is also valuable in 
giving a physical picture of the action of the control. 

When the rudder reaches the angle called for hy the control, a 
check must he made to find whether the system T^rill thereafter operate 
linearly or whether the maximum rate called for will again exceed that 
available. The rate of rudder motion called for hy the control is 
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Because the maximum rate available is E, this condition can he 
satisfied only when 


= 


E 


A pair of grid lines determined hy this relation = i - J 

is shown in figure 3- The area between these lines is designated as 
the "linear hand. " Whenever the rudder reaches the angle called for 
hy the control at a point within the area hounded hy these two lines, 
the control operates linearly, and the graphical analysis must he 
continued with the spiral for = 0.2. Whenever the trajectory 
passes outside the area hounded hy these two lines, however, the 
rudder is moving at its maximum rate and the analysis must he made hy 
using the spiral for C = 0.05- Note that the trajectory may pass 
through the linear hand without the control operating linearly. Only 
when the rudder reaches the angle called for hy the control within the 
linear hand does the contral operate linearly. 


The trajectory for the first Interval of the motion is seen from 
figure 3(a) to end outside the hand of linear operation. The rudder 
therefore reverses and moves at its maximum rate in the negative 
direction for the next inteival. The spiral is therefore shifted so 
that its origin is located at the intersection of the grid lines 
1)%^ = 0 and Dir = -EN = -( -0.125)(0.3) = 0.0375* The spiral is then 
rotated to pasa through the final point af the first interval, which 
becomes the initial point of the second interval. Values of 0 and s 
for this inteival, starting from zero at the start of the interval, 
must he determined. These values of 9 and s are used in determining 
the values of ijr for this interval and the end point of the interval. 
Finally, the values of s in the second interval are added to those 
in the first interval to determine the over-all time elapsed. The 
end of the second interval is determined in the same way as the end of 
the first interval. 


The analysis is continued in this manner for several intervals, 
the origin of the spiral being shifted each time the rudder motion 
reverses. As shown in figure 3(a) ^ the end point of the fifth interval 
occurs within the hand of linear operation. At this point, the work 
must he continued on the spiral for = 0.2. In transferring the 
work to the spiral for the linear mode of operation, allowance must, 
in general, he made for the fact that this spiral is based on 
equation (T)^ which is exp>ressed in terms of the time unit s^ (eq. (6)). 
For this reason, the derivatives D^ijr and DV and the limits of the 
linear hand must he converted to the new time unit whenever the work is 
continued on the spiral for the linear range. The required relations 
are 
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Dc^ijr 






1 + 


Stjr 


Dc^ 



Time as measured on this spiral must be converted back to units of s 


p 

when plotting the results. In the present example, s = Sg, = 


and Dc^lf = Di)r because 
has no effect in this case. 




is zero. The transformation therefore 


The spiral for =0.2 is located on its own grids so that 
its origin is at D^''^ = 0, D'<lf = 0. Values of and Dilf from the 

last point are replotted on these grids and the spiral is rotated to 
pass through this point, as shown in figure 3(b). The work then 
proceeds as before. 


In the example chosen, the motion almost immediately passes out- 
side of the band of linear operation again. The analysis is therefore 
shifted back to the spiral for ^ = 0.05 for the seventh interval. 

This interval again ends in the linear range, and thereafter the motion 
remains linear and danips out exponentially. 


Time histories of the yaw angle the control function 

and the rudder angle 5^., obtained from the preceding analysis 
plotted in figure 4. 


, are 


Example II .- The following parameters are assumed: 

36 ^ 




= 1 


? = 0.05 (oj = 0.9987) 

N = 0.3 


■•max 


= ±0.125 


from which may be derived, from equation (8), 


= 0.2 (ut = 0.9798) 
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The following initial conditions are assumed: 

^0 = -0*575 

(rftl^)o = 0 

Since the rudder is constrained to move in proportion to the yawing 
velocity as long as it is not against a stop, 

8ro =0 


Prom equation Cn) the following equation may be derived: 

= 0.575 

The graphical work for this problem is shown in figure 5- This example 
illustrates a procedure not found in example I in that the angle of yaw 
may be read fTOm suitable sets of vertical grid lines whenever the 
rudder is against a stop. These grid lines cannot be used when the 
rudder is moving at a constant rate. Initially, the control is 
operating in its linear range; therefore, the spiral for = 0.2 is 

used for the first interval (fig. 5(a)), with its origin located at 
D^ijr = 0, Ihlr = 0. The ndder always reaches its stop at a given value 

of yawing velocity, because Sj. = Ihlf — - 


dot 


The values of Dijf where the 


rudder reaches its stops are 


D* = ± ■_ = ± Q-. 125 - = ±0.125 

85r/8lh|r 1 


These lines are shown as limits of the linear band in figure 5* When 
the trajectory passes outside one of these lines, the work must be 
continued on the spiral corresponding to 5 = 0.05. This spiral is 
centered on the grid lines = 0, Dijr = 0, as shown by equations (20) 

and ( 21 ). Values of s may be determined at each point on the trajec- 
tory as described for example I. If desired, values of ijr may be 
obtained from the values of and Bf, as described for exanple I, 

or they may be read directly from STiltable sets of vertical grid lines 
constructed in accordance with equation ( 23 ). This equation shows that 
the line for ijr = -W^rQ passes throu^ the origin of the spiral. When 

the trajectory is above the line D\lr = 0.125, therefore, the grid line 
■'1^ = -( 0 . 3 ) ( 0 . 125 ) = - 0.0375 passes through the origin. When the 
trajectory is below the line Dilr = -0.125, the grid line ijr = 0.0375 
passes through the origin. When the system is operating in the linear 
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range, of coirrse, the grid line T|f = 0 passes throtigh the origin. 
The grid lines corresponding to these conditions are shown by dotted 
lines in figure 5* 


Time histories of the yaw angle \}r, the control function 

and the rudder angle Sj. obtained for this problem are plotted in 
figure 6. 


36j. 


Exan^)le III .- The following parameters are assumed: 

hb-r 

^ = 0.05 ((s = 0.9987) 

N = 0.3 
R = to. 125 


from which may he derived, from equation (8), 

= 0.0395 Co)c = 0.9992) 


The following initial conditions are assumed: 

1^0 = -0.5 

(Dt)o = 0 



- 0.25 


from which may he derived, from equation CH)^ 

(D2iir)Q = 0.575 

The graphical work for this problem is shown in figure Y, curve 1. 
The procedure required is very similar to that of exanple I, althou^ 
different equations are required to determine the end points of the 
intervals and the limits of the lineeir band. 


The end of an interval is reached when 


= 5_ + Rs 

^ ^0 



( 24 ) 


From equation (22), the value of ijr may be obtained as follows: 
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ilr = ^ - N^8r^ - 2 ^r) - ENs (25) 

Substitute equation (25) into equation (24), replace s by s/cu, and 
solve for x. The result is 


X = eR^N + A + 5 ^ cd/n + - 2^mRN (26) 

Numerical values may be substituted into equation (26) to obtain x as 
a function of 9. A curve of x as a function of 6 satisfying this 
relation nay be plotted on the graph, as shown in figure 7 hy a dashed 
line. The intersection of this curve and the spiral determines the 
time at which the rudder reaches the angle called for by the control. 
When the values of x and 9 are plotted, it must be remembered that 
the X- and y-axes are considered to pass through the origin of the 
spiral, not the origin of the D^i|f, grid as shown in the figures. 

Furthennore, values of 9 are measured about the origin of the spiral. 


The complete solution may now be carried out as it was in the 
previous examples. This solution, shc^m in figure J, c\irve 1, indicates 
that the motion approaches a steady htmting oscillation at an amplitude 
sli^tly smaller than the initial angilltude. Time histories of the 


control function "the rudder angle 5p, obtained from the 

preceding analysis, are plotted in figure 8(a). 


In order to show the behavior of the system when, the motion is 
started with smaller initial amplitudes, curves 2, 3^ and 4 of figure 7 
were plotted. The resultant time histories are shown in figures 8(b) 
to 8(d). These results indicate that, if the motion exceeds the linear 
range only slightly, eis in curve 3, the oscillation continues to build 
up ■until the steady hunting condition of curve 1 is approached. On the 
other hand, if the oscillation starts at a slightly smaller amplitude, 
as in curve 4, it remains stable and damps out. Note that in the 
solution for the linear range of operation (fig. 7(^3)) a transfonmation 
of the time scale is required. 


Example IV .- The following parameters are assiimed: 

W " ^ 

^ = 0.05 CcD = 0.9987) 

N = 0.3 

Bj. = ±0.250 
•^max 
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from which may he derived, from equation (Q), 

= 0-0395 (o^ = 0.9992) 

The following initial conditions are assumed: 

^0 " 0.575 

(Dt)o = 0 


The control attempts to move the rudder in accordance with the relation 

35 j, 

5 r = Initially, the control calls for a value of 6 ^ equal to 

(- 0 . 575 )( 2 ) = -I.15, which is greater than the maximum rudder deflection. 
The initial rudder angle is therefore the maximum value. 


6ro = -0.25 

From equation (ll) the following equation may he derived: 

=0.650 

The graphical work for this problem is shown in figure 9 - Time histories 

of the control function and the rudder angle plotted in 

figure 10. 

Example V. - The following parameters are assumed: 




36 j . 


0.05 (co = 0.9987) 


N = 0.3 

R = ± 0.2 


max 


±0.25 


from which may he derived, from equation ( 8 ), 


= 0.158 (cDc = 0.9874) 
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The following initial conditions are assumed: 


Tiro = -0.325 

(Dt}t)q = 0 


5ro = -0*25 

fix)m which may he derived, from equation (ll), 

(D2ilr)Q = O.k 


In this problem, the graphical work for which is shown in figure 11, 
curve 1, somewhat more complicated relations exist for determining 

the limits of linear operation. The quantity V + Dif determines 

the rudder position called for by the controlling mechanism. Initially, 
this quantity equals -O.65O, whereeis the rudder is against its negative 
stop at -O.25. The rudder therefore stays at its negative stop for the 
first intei*val. The rudder leaves its stop when the control function 





+ Dt = -0.25 


This relation determines a straight line on the \lr and Dijr grids 
which may be drawn in as shown in figure 11. A similar line is drawn 
to show when the rudder leaves its positive stop. 


After the rudder leaves its stop, a check must be made to determine 
whether the rudder will operate in the linear range. The rate of rudder 
motion called for by the control is 


D5j. 


Di — ^ + 1)2^ — ^ 
SDilf 


A second pair of lines determined by this relation, 

- 0,2 = Di}((2) + I)2ilf(l), is shown in figure 11. The area bounded by the 
two pairs of lines is designated as the "linear region" in the figure. 
Since the first interval ends outside the area bounded by these lines, 
the rudder moves at its maximum rate in the positive direction for the 
second inte:rval. The solution proceeds as in case I or III. 


The second Interval ends when the rudder reaches its maximum 
deflection or when it reaches the angle called for by the control, 
whichever occiors first. The time that the rudder must travel to reach 
its positive stop after the start of an interval is determined from the 
relation 
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&i» + Rs = 

u max 

from which 

c = i = 

0) R 

The value of 9 at which this condition is satisfied, in the present 
instance, is 

9 = ( 0 . 9987 ) = 2.49 radians = 14-3° 

0.2 


The conditions under which the rudder reaches the angle called for by 
the control are given by the formula 



In this formula, 9 q is the angle between the positive x-axis and 
the initial point of the interval, measured clockwise. This equation 
is an extension of equation (26) to the present problem and is derived 
in a similar manner. A few points of a curve of x as a function of 
9 satisfying this relation may be plotted on the graph. The inter- 
section of this curve and the spiral determines the time at which the 
rudder reaches the angle called for by the control. Alternatively, 
a plot of the rudder angle and the control function as functions of s 
may be made to determine this point. 

In the present Instsince, the rudder reaches its positive stop 
before it reaches the angle called for by the control. In the third 
interval, therefore, the rudder is against its positive stop. The 
solution proceeds as it did for the first interval. 


Continuation of the solution, shown in figure 11 , curve 1 , results 
in a closed cycle. The system therefore performs a steady hunting 
oscillation for the initial conditions chosen. Time histories of the 

and the rudder 

angle 5 j. obtained from the preceding analysis are plotted in 
figure 12(a). 


yaw angle ijf , the control function 


I 


+ Dilr 
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In order to show the behavior of the system when the motion is 
started with a smaller initial amplitude, curve 2 of figure 11 was 
plotted. This solution differs from that of curve 1 in that at the 
end of the first interval the rudder reaches the angle called for hy 
the control before it reaches its positive stop. The details of the 
graphical work are somewhat similar to those required for example I. 
Time histories of the motion are shown in figure 12(b). 


DISCireSION OF EESUITS 


Althou^ the preceding discussion is intended primarily to 
describe a graphical method for determining the motion of an auto- 
matically controlled system when saturation effects are present, a 
few results of interest may be obtained from the examples presented. 

Variation of stability with amplitude .- First, some general 
statements may be made regarding the vairiatlon of the stability of a 
system of the type considered as a function of amplitude. At 
amplitudes within the linear range of operation of the control, the 
period and damping may be predicted by well-known methods. For 
practical controls, of course, the dan 5 >ing should usually be considerably 
greater than that of the tincontrolled system. At very large amplitudes , 
far beyond those where the linear range is first exceeded, the period 
and damping of the system approach those of the uncontrolled system. 

This tendency results from the fact that the moments supplied by the 
control, when either its deflection or rate is limited, are very small 
compared with the inh erent restoring and damping moments acting on the 
system at sufficiently large amplitudes. At intermediate amplitudes, 
the behavior of the system depends on how the saturation effects 
Influence the phase relationship between the notion of the system and 
that of the control. If only the control deflection is limited, as in 
example II or IV, the phase relationship tends to be unaffected, and a 
smooth transition occurs between the characteristics of the controlled 
system and that of the uncontrolled system as the amplitude is increased. 
If only the control rate is limited, however, as in example I or III, 
a phase lag approaching 90 ° "the motion of the control is introduced 
and, as a result, the damping at intermediate an 5 >litudes is decreased. 
Whether actual Instability occnirs in some amplitude range depends 
primarily on the amount of damping present in the uncontrolled system 
and in the controlled system in its linear range of operation. In 
example III presented previously (fig. 8), instability occuired at an 
amplitude slightly beyond the linear range and continued up to several 
times this aa^jlitude. With both rate and deflection limited, a phase 
lag Intermediate between those produced by the rate and deflection 
limitations may be present. If the deflection range is sufficiently 
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large, however, the rate limitation may prevent the control from ever 
reaching its stops; therefore the deflection limitation becomes 
unimportant . 

Examples I and II, presented previously, in which the control 
produced primarily an increase in damping, show that the control 
remained effective in producing damping even at amplitudes several 
times that at which the saturation effects were first encountered 
(figs. 4 and 6). This result may he of importance in the design of 
yaw dampers for airplanes because it indicates that, if some decrease 
in danping at the larger amplitudes can be tolerated, the power output 
of the controlling device may be made considerably less than that 
required to provide linear operation throughout the range of amplitudes 
encountered. As was previously mentioned, flight tests have shown that 
increased damping of the lateral oscillations of aiiplanes is required 
mainly at small amplitudes. 


Extension of the examples to other conditions .- Some reduction in 
the number of parameters in the problem has already been made by express- 
ing time in terms of the nondimensional parameter s . The results of 
any given problem may therefore be applied to systems with any natural 
frequency by converting back from the nondimensional to the real time 
unit. A further reduction in the ntunber of parameters is possible by 
dividing throTigh the equations for the variations of t and Di|r 
(eqs. (9) and (lO)) by the quantity EN. The problem is then completely 

specified in terms of the qioantities C, &r/^^ ^/M, Df/BN, N 

S6j, 

and N A given solution may therefore be applied to systems with 


other values of R and N provided the initial conditions are changed 
so that the quantities b^Q^R, ^q^RN, and Dil<QyRN remain the same and 

provided the gearing constants are changed so that the quantities R 
and N remain the same. Values of Sr/R, ''1''/RN, and Eft/RW 


obtained from the given solution apply to the new system. The values 
of bp, ilr, and D'>lf may then be obtained by multiplying these parameters 
by the appropriate factors. 


Application to actual automatic controls .- Although actual auto- 
matic controls do not follow the controlling quantities with zero lag 
as assumed in the foregoing analysis, a method is suggested whereby 
the same type of analysis may be applied to give approximate solutions 
for saturation effects on actiial automatic controls. For this puip)ose, 
the gearing constants 9br/9''lf and dbr/SC'l^ should be selected to give 
the same amplitude ratio and phaise angle between the control motion and 
the motion of the system as exists on the actual automatic control at 
the frequency at which the controlled system oscillates in its linear 
range. Usually this frequency may be approximately determined from 
elementary considerations. The applicability of this method depends 
on the condition that the satiiration effects do not change the frequency 
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very mucli. This condition may he expected to exist when the restoring 
moment on the uncontrolled system is fairly large compared with the 
restoring moment supplied hy the control. The accuracy of the results 
obtained hy this method prohahly depends to some extent on the actual 
details of the mechanism under consideration. This method is limited 
in applicability to cases in which the controlled member can be 
approximated as a single-degree-of-freedom system. No limitation is 
placed on the complexity of the controlling mechanism itself , inasmuch 
as experimentally determined frequency-response data could be used to 
select the gearing constants. 


CONCLUSIONS 


A graphical method is presented for determin in g the motion of a 
freely oscillating system of one degree of freedom stabilized by a 
controlling device which applies control force in proportion to the 
displacement of the system, to its rate of change of displacement , 
or both. The controlling member is asstnned to have limitations on 
its maximum deflection, on its maximum rate of movement, or both. 
Although the report is intended primarily to describe this -graphical 
method, the following concltisions may be obtained from the examples 
presented: 

1. At sufficiently small amplitudes the period and damping of 
the system correspond to those provided by linear operation of the 
control whereas, at very large amplitudes, the period and damping 
approach those of the uncontrolled system. At intermediate amplitudes, 
if the deflection of the control is limiited, a smooth transition 
between these two conditions of period and daumping takes place whereas , 
if the rate of control movement is limited, the damping may be reduced 
below that of the uncontrolled system. In some cases, limiting the 
rate of control movement may produce Instability over a range of 
amplitudes . 

2. If the control produces primarily an increase in damping, 
the control remains effective in producing damping even at amplitudes 
several times that at which saturation effects are first encountered. 
This effect may be useful in reducing the power reqiilrements of yaw 
dampers for airplanes . 


Langley Aeronautical Laboratory, 

National Advisory Commiittee for Aeronautics, 
Langley Field, Va., August l8, 1953* 
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